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Abstract 

It is shown that, by letting wavenumbers and frequencies complex in Hirota's bi- 
linear method, new classes of exact solutions of soliton equations can be obtained 
systematically. They include not only singular or A^-homoclinic solutions but also 
A^-wavepacket solutions. 
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Recent awareness of homoclinic solutions of soliton equations [1,2] has made it 
possible to apply them for studies of chaotic phenomena in nonlinear, infinite 
dynamical and complex systems, e.g., [3]. In contrast to the low-dimensional 
ordinary differential equations, partial differential equations may have multiple 
or infinite homoclinic structures. Therefore, it is important to find and classify 
homoclinic solutions systematically. 

The key idea of this paper is to make wavenumbers complex in soliton so- 
lutions. Usually this procedure fails because the dependent variable becomes 
complex for real- valued equations. However, if we keep the complex conjugate 
pair for wavenumbers, the difficulty can be avoided. 

Let us start with the simplest system, the KdV equation, 

Ut + 6mM^ + Ua:xx = 0, (1) 



where Ux = d^u = du/dx. Letting u{x,t) = 2(\n f{x,t))xx, the bilinear form 
becomes [4] 

M/ ■ / = 0, (2) 
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where 



M = M{Dt,D,) = D,{Dt + Dl), 



(3) 



and is the Hirota's bihnear operator, 

Dxf -g^id:^- da:')f{x)g{x%'=^. (4) 

The iV-sohton solution is given by 

(N N \ 

j=i j>i I 

with 

r]j = kjX + ^jt + ?7oj, (6) 
Q,- = (7) 



where I]^=o,i indicates the summation over all possible combinations of //j = 
0,1 for i = i,---,Ar. 

This solution gives solitons for real kj. If we make complex conjugate as 
kj-x — k* for j = 2, 4, ■ ■ ■ , 2M and real for j = 2M + 1, - ■ ■ ,N, it no longer 
gives classical solitons. We note that the proof of the A^-soliton solutions by 
mathematical induction [4] is still valid for complex wavenumbers. 

For M = 1 and L = N — 2M — 0, the function / and the solution can be 
described by 

/ = l + 2exp^iCosCi + ai2exp(2^i), (9) 

, r^2^^^^ cos(Ci + 2(pi) + ai2C^i cos(Ci - 2(pi) + 2[ai2 cos^ (pi - sin^ ipi] 
u — 4A^ 



(e-€i + ai2e«i + 2cosCi)2 



where and Q are the real and imaginary parts of rji, and 



Fig. 1. The 3D plot of a one-singular-point solution of the KdV equation with 
ki = exp(i7r/4). The constants ijio are zero throughout the illustrative examples in 
this paper. 

for ki = Kiexpiipi. Since is negative, it can be shown that this solution 
possesses one singular point (/ = 0) in a; for given t. Figure 1 shows the 3D 
plot of this solution. 

It looks as the singular point travels at the velocity —Qij./kir (subscript r 
denotes the real part) with oscillating back and forth in the stripe surrounded 
by lines = [ln(v^l — 012) ± !]/(— 012), which are derived by / = and 
cosCi = ±1 in (9). The singularity behaves as u Ug — —2{x — Xq)"^, which 
is derived by a simple zero of / in x. It should be noted that Ug is an exact 
steady solution of the KdV equation. It is also observed that a small pulse with 
the sign opposite to the singularity is travelling with the velocity —flu/ku. 

The interaction between two singular points for (M, L) = (2, 0) is illustrated 
in Figure 2. The asymptotic behavior of the singular point for ki can be given 
as follows. In the moving frame Xi — x + {flir/kir)t, rj^ can be written as 

r)3 = kzrXi + a^rt + iikziXi + a^it) + rj^o, (12) 
where — f^3r — k^r^ir/kir and a^i — fl^i — ku^ir/kir- H crsr < 0, 



/ ~ 1 + e^i + e^2 + 0126^1+"^ = 0, 



(13) 
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Fig. 2. The 3D plot of the interaction between two singular points of the KdV 
equation with ki = exp(i7r/4) and k2 = 2^/^ exp(i7r/6). 

as i — > oo. Similarly, in the moving frame — x + {^sr/ksr)^, 

r]i = kirX^ + airt + iikiiXs + ant) + rjio, (14) 

where air — —kircr^r/ksr and an — Qij — kuQ^^r/k^r- Assuming kir/ksr > 0, 
we have air > and 

/ ~ aise^^+^^il + 6"^+^' + e''*+'^* + 0346^'^+''*+^^+^*], (15) 

as t ^ 00, where the complex phase shift is given by 63 — 51 — In(ai3a23). 

In the hmit t —00 and with the same assumption as the above, the moving 
frames Xi and X3 lead to the asymptotic forms (15) and (13) with replacement 
(1, 2) ^ (3, 4) and the phase shift 51 = ^2 = 111(013014). From the definition of 
ttij, we have 5i — 5s and the conservation law of the total phase shift holds, 
i.e., 5i + 52 — 5s + 5^. The conservation law does not depend on the signs of 
kir/ksr and air- 

Figure 3 shows the interaction between a singular point and a soliton for 
(M, L) = (1, 1). Under the same assumptions of the signs of kir/k^r and air, 
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Fig. 3. The 3D plot of the interaction between a singular point and a soliton of the 
KdV equation with ki = exp(z7r/4) and ^2 = 2^/^. 

the asymptotic behavior of the sohton is 



^3r„„ 1 2^ 



u - ^sech"^, (16) 
2 2' ^ ' 



as t 



-oo and 



u ~ 



%sech^ 



+ 53 



(17) 



as t — s> cxD, where the phase shift is given by ^3 = In(ai3a23)- Similarly, the 
behavior of the singular point is determined by the zero of (13) and (15) with 
(1,2) <-> (3,4), the factor in front of the square bracket in (15) replaced by 
e''^, and the phase shift 51 = ^2= Inaia as i — > 00 and t — > —00, respectively. 
Again, the total phase shift is conserved. 

As we see in the above examples, the singular points behave as if they keep 
their identities after they collide with each other like solitons. Therefore, they 
may be called singulons. 

At first sight, it seems that this singularity makes the solutions physically 
meaningless. However, the singularity shown above does not always appear. 
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As such an example, the Boussinesq equation of the form 



(18) 



is considered. Substituting u = 2{]iif)xx, the bihnear equation [5] is given by 
(2) with the operator 

M^D^-Dl- Dl. (19) 

The solution is given by (5) with the dispersion relation 

VLi = ±ki^Jl + kl (20) 

We note that this solution with complex wavenumbers corresponds to a certain 
limit of the periodic solutions by two-dimensional ^-function (3.1) of [6] with 
r]i = r)2, although it is implicit that the wavenumber and the frequency may 
be complex in [6] . 

If wc choose pure imaginary wavenumbers kj^i = —kj = iKj^i and Kj > 1 for 
j = 2, ■ ■ ■ , 2M and L = 0, we obtain so-called homoclinic solutions, analogous 
to homoclinic orbits in ordinary differential equations. Strictly speaking, they 
are temporally homoclinic and spatially periodic (or quasi-periodic) solutions. 

As an example, let us consider the simplest complexation, i.e. (M, L) = (1, 0). 
The asymptotic behaviors of the solution given by (10) are 

4KlexpCiCOs{Ci + 2ifi), (21) 

as t ^ — oo and 

u ~ 4iriVai2exp(-ei) cos(Ci - 2ip^), (22) 
as t — > oo, assuming that flir is positive. 

If we put ki — —k2 — iKi and Ki > 1, the frequency becomes real and can 
be chosen as Qi = Q2 = Ki \J — 1 . In this case, the interaction coefficient 
is ai2 = (4i^i — l)/{Kf — 1) and positive. Since ai2 > 4 and / can be written 

as 

/ = 2V^cosh(^i + (5i) + 2cosCi, (23) 
with 5i — (In 012)72, this solution remains finite for all x and t. 
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Fig. 4. The 3D plot of the one-homoclinic solution of the Boussinesq equation with 
ki = 2i . 

As t — >■ —oOj the solution tends to a homoclinic point u = and behaves 
as w ~ — exp ^1 cos Ci, where = flit + 3?(?7io) and (i = Kix + 3(?7io). 
Similarly, as t — > oo, ~ —AKf/ai2 exp(— ^i) cos Ci- This behavior is no more 
than the linear approximation of the Boussinesq equation around u = 0. Here 
ai2 plays a role of shift of time and the offset of the amplitude. Figure 4 shows 
the 3D plot of the one-homoclinic solution. 

One remarkable feature of the Boussinesq equation with the spatial period 

Ix comparing with the nonlinear Schrodinger equation [1] is that there is a 
countably infinite number of homoclinic orbits due to the short wavelength 
instability kj = iKj with Kj = 2t{ j/l^ > 1. 

The necessary and sufficient condition for the regularity of the {M,L) — (1,0) 
solution is ai2 > 1. Using 



ai2 = 



1 + K^e^^^l - 1 - AKl cos 2^1 + SXf 
1 + Xfe^^^l - 1 - AKl cos 2(^1 - 3Xf ' 



(24) 



and > 1, the condition can be summarized as 



3 1 



(25) 



cos2(/?i < — - — 



4 AKl ' 



For a complex wavenumber ki, this solution may no longer be called homo- 
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Fig. 5. The 3D plot of the one-wavepacket solution of the Boussinesq equation with 
a complex wavcnumbcr ki = 2 exp(i(/?i), (y^i = cos^-'^(l/4) . 

clinic. Figure 5 shows the 3D plot of the solution with a complex wavenumber 
ki satisfying the condition (25). 

The solution shows that a wavepacket travelling at the group velocity —Vtir/ kir 
of waves with the phase velocity —Qir/ku- Therefore, it may be called a one- 
wavepacket solution. Similarly, we may construct two- and A^-wavepacket so- 
lutions, and interactions between solitons, homoclinic solutions and wavepack- 
ets. 

We have shown that the Boussinesq equation has richer geometrical structures 
in phase space than the KdV equation, although the homoclinic solutions may 
be far from modelling shallow water waves because of the assumption Ut ~ —u^ 
in the derivation of the equations [7]. The complexation of wavenumbers in 
solitons shown in this paper has a wide application and is expected to give a 
tool for finding a new class of exact solutions. 

References 

[1] B. M. Herbst and M. J. Ablowitz, Phys. Rev. Lett. 62 (1989) 2065. 
[2] M. Umeki, RIMS kokyuroku 974 (1996) 

[3] A. R. Bishop, M. G. Forest, D. W. McLaughlin and E. A. Oberman, Phys. Lett 
A 126 (1988) 335. 



8 



[4] R. Hirota, Phys. Rev. Lett. 27 (1971) 1192. 
[5] R. Hirota, J. Math. Phys. 14 (1973) 810. 

[6] A. Nakamura, J. Phys. Soc. Jpn. 47 (1979) 1701. 
[7] J. W. Miles Ann. Rev. Fluid Mech. 12 (1980) 11. 



9 



